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■ We report the first observation of e + e — annihilations into hadronic states of 
' positive C-parity, p°p° and </>p°. The angular distributions support two-virtual- 
i photon annihilation production. We also report the observations of e + e _ — » <f>r) 

and a preliminary result on e+e~ — > p+p — . 
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, 1. Introduction 

m ■ 

■ The large datasets collected by the B factories provide unique opportu- 
, nities for studying rare processes and discovering new states. We report 

several observations of e + e~ annihilations into quasi-two-body hadronic fi- 
C ■ . nal states with C — ±1 at BABAR [1-3]. A new avenue for the study of 

hadron production mechanisms is opened with these observations, and a 
testing ground for QCD at the amplitude level is provided. 



X 



2. ev-pv.y 



The process e + e~ — * hadrons at center-of-mass (cm.) energy i/s far below 
the Z° mass is dominated by annihilation via a single virtual photon, thus 
yielding final state charge-conjugation parity C = — 1. The Two- Virtual- 
Photon- Annihilation (TVPA) process, depicted in Fig.[lJ with positive final 
state C parity, has been ignored in incorporating the total hadronic cross 
section in e + e~ annihilations into calculations [4] of muon g— 2, and the 
running of the QED coupling constant, a. 

The present analysis uses a 205 fb _1 data sample collected at the T(4S f ) 
resonance, and 20 fb _1 collected at cm. energy 40 MeV lower, using the 
BABAR detector at the SLAC PEP-II asymmetric-energy e + e" collider. 
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The BABAR detector is described in detail elsewhere [5]. 

Events with four well-reconstructed charged tracks and net charge zero 
are selected. The \ 2 probability of the fitted four track vertex is required 
to exceed 0.1%, and two oppositely charged tracks must be identified as 
pions; the other pair must be identified as two pions or two kaons. We 
accept events with four-particle invariant mass within 170 MeV/c 2 of the 
nominal cm. energy. Loose signal regions are defined by the mass ranges 
0.5 < m w+7r - < 1.1 GeV/c 2 and 1.008 < m K+K - < 1.035 GeV/c 2 . The 
extracted p°p° and (f>p° yields in these intervals are 1243 ± 43 and 147 ± 13 
events, respectively. 

The efficiency-corrected production angular distributions are shown in 
Fig. [2] where 9* is defined as the angle between the p° (4>) direction and 
the e~ beam direction in the cm. frame. The observed sharply peaking 
| cos 9* | distributions are consistent with the TVPA expectation [6] , which 
is approximated by: 

da l + cos 2 #* 
d cos 9* * 1 - cos 2 9* ' ( ' 

For the signal mass regions defined above, and | cos#*| <0.8, we obtain 
the following results for the TVPA cross sections near y/s = 10.58 GeV: 

cr fid (e + e~ -> p°p°) = 20.7 ± 0.7(stat) ± 2.7(syst) fb 
CTfi d (e + e- -> (j)p a ) = 5.7 ± 0.5(stat) ± 0.8(syst) fb. 

The measured cross sections are in good agreement with the calculations [6, 
7] . The Standard Model calculations of the anomalous magnetic moment of 
the muon and of the QED coupling constant rely on measurements of low- 
energy e + e _ hadronic cross sections, which are assumed to be entirely due 
to single-photon exchange. We have estimated the effect due to the TVPA 
processes [6] and find it to be small compared with the current precision [4] . 




Fig. 1. The two- virtual-photon annihilation diagram. 
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3. e+e — > (fir) 

The process e + e~ — > J/iprj c and other double charmonium processes are 
observed [8] at rates approximately ten times larger than expected from 
QCD-based models [9]. Various theoretical efforts to understand this have 
been made recently [10] . An alternate avenue of investigation is provided by 
the process e + e~ — > tfirj, which also involves a vector — pseudoscalar (VP) 
final state. Different models predict different s dependences for the cross 
section, and so it is interesting to investigate this by comparing a mea- 
surement at y/s — 10.58 GeV to the CLEO measurement at y/s = 3.67 
GeV [11]. 

This analysis uses 204 fb _1 of e + e~ colliding beam data collected on the 
T(45) resonance at y/s = 10.58 GeV and 20 fb" 1 collected 40 MeV below. 
Events with exactly two well-reconstructed, oppositely charged kaon tracks 
and at least two well-identified photons are selected. We fit the two tracks 
to a common vertex, and require the \ 2 probability to exceed 0.1%. Each 
photon candidate is required to have a minimum laboratory energy of 500 
MeV. Events with a reconstructed K + K~-f~f invariant mass within 230 
MeV/c 2 of the e + e~ cm. energy are accepted for further study. 

We define the mass window as 1.008 < tukk < 1.035 GeV/c 2 , and 
extract 24 ± 5 cfrr) signal events in the <fi mass window, with r\ — > 77. The 
significance is estimated to be 6.5 sigma. 

The final radiation-corrected cross section for 1.008 < < 1.035 

GeV/c 2 within | cos f5* | < 0.8 near ^ = 10.58 GeV is: 

cr fid (e + e~ (jyq) = 2.1 ± 0.4(stat) ± O.l(syst) fb. 
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The cross section, extended to the full range of cos 9* by assuming a 1 + 
cos 2 8* distribution, is: 

cr(e + e" (jyrj) = 2.9 ± 0.5(stat) ± O.l(syst) fb. 

There is currently no direct prediction for the cross section of this process 
at this energy, but the e + e~ — > VP cross section is expected to have a 
1/s 2 [12] or 1/s 4 [13,14] dependence in QCD-based models. A comparison 
between our result and that of CLEO, (rj = 2.1±{% ± 0.2 pb) at ^/s = 3.67 
GeV (continuum) [11], favors a 1/s 3 dependence (Fig. [3]). 
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Fig. 3. Cross section extrapolations based on BABAR's measurement at ^/s = 
10.58 GeV assuming a 1/s 3 (black) or 1/s 4 (red) energy dependence. The bands 
show one standard deviation uncertainties in the extrapolations. The CLEO mea- 
surement at ,/s = 3.67 GeV is also shown. 



4. e+e — > P~*~P (preliminary result) 

Since charged p's are involved, the e + e~ — * p + p~ process is unlikely to 
occur through TVPA [1,6,7], unless there is significant final quark recombi- 
nation between the products of the two virtual photons, or unless there is 
significant final state interaction (e + e~ — > p°p° — > p + p~ ) [15]. Assuming a 
one-photon production mechanism, this VV (p + p~) final state is described 
by three helicity amplitudes. A study of this reaction can then provide an 
experimental test of QCD at the amplitude level [14] [13] through investi- 
gation of the final states angular correlations. 

This analysis uses 343 fb -1 of e + e~ colliding beam data collected on the 
T(45) resonance at *Js = 10.58 GeV and 36 ftr 1 collected 40 MeV lower. 
Events with exactly two well-reconstructed, oppositely charged tracks iden- 
tified as pions and at least two well-reconstructed 7r°s are selected. We fit 
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Fig. 4. The s-weighted and efficiency corrected a) cos f?+ b) cos 9— c) tp+ d) </?_ e) cos 6* 
distributions for e+e — — > p+ p~ . The magenta dashed curves show the contributions from 
Foo, the blue dotted curves are Fin, the blue dashed-dotted curves are Fu, and the solid 
red curves show the total result. 



the charged tracks to a common vertex, and require the y 2 probability to 
exceed 0.1%. Each ir° is reconstructed through its 77 decay channel by 
requiring the two photon invariant mass to be within the range [0.1,0.16] 
GeV/c 2 , and then constraining its mass to the nominal value. We accept 
events with \m„+„o„- n o — E cm \ < 0.28 GeV and \Ap\ < 0.2 GeV/c, where 
E cm is the total cm. energy, and Ap is the momentum difference between 
the 7r + 7r°7r~7r° system and the e + e~ system. 

We define the p mass intervals as [0.5,1.1] GeV/c 2 , and extract 308±25 
p + p~ signal events in the defined mass region. The significance is estimated 
to be 9.5 sigma. 
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Assuming p + p is produced through one photon or T(45), there are 
three independent hclicity amplitudes (F^,,, pjv is the helicity of p + / p~), 
F 00 , Fio, and F n (F 10 = F_ 10 = F 0±1 , F u = F_ x _i) [16]. The one- 
dimensional projections for the decay angles involved can be expressed as: 

f '' A oc (sin 2 r|F o| 2 + 2(l + cos 2 r)|F 10 | 2 + 2sin 2 r|F 11 | 2 ) (2) 



dcos9* 
dN 



d cos 0± 

dN 



cx (cos 2 c?±|F 00 | 2 + (1 + cos 2 0±)|Fi„| 2 + sin 2 9±\F n \ 2 ) (3) 
oc (|F 00 | 2 + (4-cos2( /3± )|Fio| 2 + 2|Fii| 2 ) (4) 



dip± 

where 9* is the p production angle, 9± (<p±) is the helicity (azimuthal) angle 
of the pion from p decay. From the two dimensional mass fit (ir + ir and 
7r~7r°), we can calculate a p + p~ signal sWeight [17] for each event (including 
those events outside the defined p mass window) and use it to produce 
signal angular distributions. We fit the five angular distributions to Equs.[21 
[3] and [4] simultaneously by minimizing x 2 ■ The correlations among the five 
angles are neglected; this is justified by means of fits to events generated 
according the assumed PDFs (toy MC). We normalize the amplitudes such 
that |.Foo| 2 + 4|Fi | 2 + 2|^n| 2 = 1 since we have 1 F 00 , 4 Fi and 2 Fn 
amplitude contributions. The normalized amplitudes from the fit are found 
to be in the ratio: |F 00 | 2 : |Fi | 2 : |Fn| 2 = 0.51 ± 0.14(stat) ± 0.02(syst) : 
0.10 ± 0.04(stat) ± O.Ol(syst) : 0.04 ± 0.03(stat) ± O.OO(syst), and |F 00 | 2 
deviates from 1 with significance more than 3 sigma. This disagrees with a 
QCD prediction [14], and suggests that either the decay is not dominated 
by single-photon exchange as naively expected, or that the QCD prediction 
does not apply to data in our energy region. The final radiation-corrected 
cross section for 0.5 < m p ± < 1.1 GeV/c 2 , and within | cos 6 1 * | < 
0.8, | cos 6*± | < 0.85, at near y/s = 10.58 GeV (assuming only one-photon 
production) is: 

CTfid(e + e~ -> p + p~) = 8.5 ± 0.7(stat) ± 1.5(syst) fb. 

We extend the cross section calculation from our acceptance region to the 
full phase space using the fitted amplitude values, and find 20.0±1.6(stat)± 
3.6(syst) ± 1.7(ampl) fb; the third uncertainty is due to the amplitude un- 
certainties. 
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5. Conclusion 

We report the first observation of e + e~ annihilations into hadronic states 
of positive C-parity, p°p° and 4>p°. We also report the observation of the 
process e + e~ — > 4>rj, and obtain preliminary results on e + e~ — > p + p~ . The 
measured helicity amplitude magnitudes from e + e~ — > p + p~ contradict a 
QCD prediction at a significance of more than 3 sigma. 
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